Genomics reveals the role of admixture in the evolution of structure among sperm whale populations within the Mediterranean Sea.

In oceanic ecosystems, the nature of barriers to gene flow and the processes by which populations may become isolated are different from the terrestrial environment, and less well understood. In this study we investigate a highly mobile species (the sperm whale, Physeter macrocephalus) that is genetically differentiated between an open North Atlantic population and the populations in the Mediterranean Sea. We apply high-resolution single nucleotide polymorphism (SNP) analysis to study the nature of barriers to gene flow in this system, assessing the putative boundary into the Mediterranean (Strait of Gibraltar and Alboran Sea region), and including novel analyses on structuring among sperm whale populations within the Mediterranean basin. Our data support a recent founding of the Mediterranean population, around the time of the last glacial maximum, and show concerted historical demographic profiles in both the Atlantic and the Mediterranean. In each region there is evidence for a population decline around the time of the founder event. The largest decline was seen within the Mediterranean Sea where effective population size is substantially lower (especially in the eastern basin). While differentiation is strongest at the Atlantic/Mediterranean boundary, there is also weaker but significant differentiation between the eastern and western basins of the Mediterranean Sea. We propose, however, that the mechanisms are different. While post-founding gene flow was reduced between the Mediterranean and Atlantic populations, within the Mediterranean an important factor differentiating the basins is probably a greater degree of admixture between the western basin and the North Atlantic and some level of isolation between the western and eastern Mediterranean basins. Subdivision within the Mediterranean Sea exacerbates conservation concerns and will require consideration of what distinct impacts may affect populations in the two basins.

Genomewide analysis of sperm whale E2 ubiquitin conjugating enzyme genes.

Marine mammals are exposed to the oxidative stress induced by hypoxia/reoxygenation cycles yet resist cellular damage. The availability of high-quality genomes promises to provide insights on how this is achieved. In this study, we considered the ubiquitinconjugating enzymes (E2) gene family, UBE2 genes, which encodes enzymes with critical roles in cellular physiology, including the oxidative stress response. The sperm whale was the first marine mammal with a chromosome-level genome, allowing the study of gene family repertories, phylogenetic relationships, chromosome gene organization, and other evolutionary patterns on a genomewide basis. Here, 39 UBE2 genes (similar to human, including 32 intact genes, one partial gene, six pseudogenes) were identified in sperm whale genome. These genes were found on 17 chromosomes and were assigned into 23 subfamilies, 16 subgroups, and four classes based on structural characteristics and functions, phylogeny and conserved domains, respectively. Although the gene structure and motif distribution of sperm whale UBE2 genes are conserved in each subfamily, motif variation and intron gain/loss may contribute to functional divergence. Segmental duplications were detected in six gene pairs, which could drive UBE2 gene innovation in the sperm whale. Contrasting seven cetaceans and five terrestrial taxa, we found that cetaceans have experienced shifts in selective constraint on UBE2 genes, which may contribute to oxidative stress tolerance during the adaptation to aquatic life. Our results provide the first comprehensive survey of cetacean UBE2 genes.

Integrated Full-Length Transcriptome and RNA-Seq to Identify Immune System Genes from the Skin of Sperm Whale (Physeter macrocephalus).

Cetaceans are a group of secondary aquatic mammals whose ancestors returned to the ocean from land, and during evolution, their immune systems adapted to the aquatic environment. Their skin, as the primary barrier to environmental pathogens, supposedly evolved to adapt to a new living environment. However, the immune system in the skin of cetaceans and the associated molecular mechanisms are still largely unknown. To better understand the immune system, we extracted RNA from the sperm whale's (Physeter macrocephalus) skin and performed PacBio full-length sequencing and RNA-seq sequencing. We obtained a total of 96,350 full-length transcripts with an average length of 1705 bp and detected 5150 genes that were associated with 21 immune-related pathways by gene annotation enrichment analysis. Moreover, we found 89 encoding genes corresponding to 33 proteins were annotated in the NOD-like receptor (NLR)-signaling pathway, including NOD1, NOD2, RIP2, and NF-kB genes, which were discussed in detail and predicted to play essential roles in the immune system of the sperm whale. Furthermore, NOD1 was highly conservative during evolution by the sequence comparison and phylogenetic tree. These results provide new information about the immune system in the skin of cetaceans, as well as the evolution of immune-related genes.

First evidence of terrestrial ambrein formation in human adipocere.

To date, the only known occurrence of ambrein, an important perfumery organic molecule, is in coproliths found in about one in a hundred sperm whales. Jetsam ambergris coproliths from the whale are also found occasionally on beaches worldwide. Here we report on the surprising occurrence of ambrein in human adipocere. Adipocere is a waxy substance formed post-mortem during incomplete anaerobic decomposition of soft tissues. Adipocere samples obtained from grave exhumations were analysed using gas chromatography-mass spectrometry (GC-MS). In addition to the typical fatty acids of adipocere, lesser amounts of ambrein were identified in the samples, in abundances similar to those of the major accompanying faecal steroids. The distribution of these compounds suggests that ambrein was produced post-mortem during the microbial decomposition of faecal residues and tissues. It is assumed that the adipocere matrix of saturated fatty acidsaided the preservation of ambrein over extended periods of time, because adipocere is stable against degradation. The association of ambrein formation in ageing faecal material, under moist, oxygen-depleted conditions, now requires more attention in studies of other mammalian and geological samples. Indeed, ambrein and its transformation products may be useful novel chemical indicators of aged faecal matter and decomposed bodies.

The first chromosome-level genome for a marine mammal as a resource to study ecology and evolution.

Marine mammals are important models for studying convergent evolution and aquatic adaption, and thus reference genomes of marine mammals can provide evolutionary insights. Here, we present the first chromosome-level marine mammal genome assembly based on the data generated by the BGISEQ-500 platform, for a stranded female sperm whale (Physeter macrocephalus). Using this reference genome, we performed chromosome evolution analysis of the sperm whale, including constructing ancestral chromosomes, identifying chromosome rearrangement events and comparing with cattle chromosomes, which provides a resource for exploring marine mammal adaptation and speciation. We detected a high proportion of long interspersed nuclear elements and expanded gene families, and contraction of major histocompatibility complex region genes which were specific to sperm whale. Using comparisons with sheep and cattle, we analysed positively selected genes to identify gene pathways that may be related to adaptation to the marine environment. Further, we identified possible convergent evolution in aquatic mammals by testing for positively selected genes across three orders of marine mammals. In addition, we used publicly available resequencing data to confirm a rapid decline in global population size in the Pliocene to Pleistocene transition. This study sheds light on the chromosome evolution and genetic mechanisms underpinning sperm whale adaptations, providing valuable resources for future comparative genomics.

Inter-individual differences in contamination profiles as tracer of social group association in stranded sperm whales.

Ecological and physiological factors lead to different contamination patterns in individual marine mammals. The objective of the present study was to assess whether variations in contamination profiles are indicative of social structures of young male sperm whales as they might reflect a variation in feeding preferences and/or in utilized feeding grounds. We used a total of 61 variables associated with organic compounds and trace element concentrations measured in muscle, liver, kidney and blubber gained from 24 sperm whales that stranded in the North Sea in January and February 2016. Combining contaminant and genetic data, there is evidence for at least two cohorts with different origin among these stranded sperm whales; one from the Canary Island region and one from the northern part of the Atlantic. While genetic data unravel relatedness and kinship, contamination data integrate over areas, where animals occured during their lifetime. Especially in long-lived animals with a large migratory potential, as sperm whales, contamination data may carry highly relevant information about aggregation through time and space.

Demography or selection on linked cultural traits or genes? Investigating the driver of low mtDNA diversity in the sperm whale using complementary mitochondrial and nuclear genome analyses.

Mitochondrial DNA has been heavily utilized in phylogeography studies for several decades. However, underlying patterns of demography and phylogeography may be misrepresented due to coalescence stochasticity, selection, variation in mutation rates and cultural hitchhiking (linkage of genetic variation to culturally-transmitted traits affecting fitness). Cultural hitchhiking has been suggested as an explanation for low genetic diversity in species with strong social structures, counteracting even high mobility, abundance and limited barriers to dispersal. One such species is the sperm whale, which shows very limited phylogeographic structure and low mtDNA diversity despite a worldwide distribution and large population. Here, we use analyses of 175 globally distributed mitogenomes and three nuclear genomes to evaluate hypotheses of a population bottleneck/expansion vs. a selective sweep due to cultural hitchhiking or selection on mtDNA as the mechanism contributing to low worldwide mitochondrial diversity in sperm whales. In contrast to mtDNA control region (CR) data, mitogenome haplotypes are largely ocean-specific, with only one of 80 shared between the Atlantic and Pacific. Demographic analyses of nuclear genomes suggest low mtDNA diversity is consistent with a global reduction in population size that ended approximately 125,000 years ago, correlated with the Eemian interglacial. Phylogeographic analysis suggests that extant sperm whales descend from maternal lineages endemic to the Pacific during the period of reduced abundance and have subsequently colonized the Atlantic several times. Results highlight the apparent impact of past climate change, and suggest selection and hitchhiking are not the sole processes responsible for low mtDNA diversity in this highly social species.

The Novel Evolution of the Sperm Whale Genome.

The sperm whale, made famous by Moby Dick, is one of the most fascinating of all ocean-dwelling species given their unique life history, novel physiological adaptations to hunting squid at extreme ocean depths, and their position as one of the earliest branching toothed whales (Odontoceti). We assembled the sperm whale (Physeter macrocephalus) genome and resequenced individuals from multiple ocean basins to identify new candidate genes for adaptation to an aquatic environment and infer demographic history. Genes crucial for skin integrity appeared to be particularly important in both the sperm whale and other cetaceans. We also find sperm whales experienced a steep population decline during the early Pleistocene epoch. These genomic data add new comparative insight into the evolution of whales.

Phylogenetic profiling and gene expression studies implicate a primary role of PSORS1C2 in terminal differentiation of keratinocytes.

PSORS1C2 is a gene located between coiled-coil alpha-helical rod protein 1 (CCHCR1) and corneodesmosin (CDSN) within the psoriasis susceptibility locus 1 (PSORS1). Here, we performed a comparative genomics analysis of the as-yet incompletely characterized PSORS1C2 gene and determined its expression pattern in human tissues. In contrast to CCHCR1, which is common to all vertebrates investigated, PSORS1C2 and CDSN are present exclusively in mammals, indicating that the latter genes have originated after the evolutionary divergence of mammals and reptiles. CDSN is conserved in aquatic mammals, whereas PSORS1C2 orthologs contain gene-inactivating frame shift mutations in whales and dolphins, in which the epidermal differentiation programme has degenerated. Reverse-transcription PCR screening demonstrated that, in human tissues, PSORS1C2 is expressed principally in the epidermis and weakly in the thymus. PSORS1C2 mRNA was strongly upregulated during terminal differentiation of human keratinocytes in vitro. Immunohistochemistry revealed exclusive expression of PSORS1C2 in the granular layer of the epidermis and in cornifying epithelial cells of Hassall's corpuscles of the thymus. In summary, our results identify PSORS1C2 as a keratinocyte cornification-associated protein that has originated in evolutionarily basal mammals and has undergone gene inactivation in association with the loss of the skin barrier function in aquatic mammals.

What influences the worldwide genetic structure of sperm whales (Physeter macrocephalus)?

The interplay of natural selection and genetic drift, influenced by geographic isolation, mating systems and population size, determines patterns of genetic diversity within species. The sperm whale provides an interesting example of a long-lived species with few geographic barriers to dispersal. Worldwide mtDNA diversity is relatively low, but highly structured among geographic regions and social groups, attributed to female philopatry. However, it is unclear whether this female philopatry is due to geographic regions or social groups, or how this might vary on a worldwide scale. To answer these questions, we combined mtDNA information for 1091 previously published samples with 542 newly obtained DNA profiles (394-bp mtDNA, sex, 13 microsatellites) including the previously unsampled Indian Ocean, and social group information for 541 individuals. We found low mtDNA diversity (π = 0.430%) reflecting an expansion event <80 000 years bp, but strong differentiation by ocean, among regions within some oceans, and among social groups. In comparison, microsatellite differentiation was low at all levels, presumably due to male-mediated gene flow. A hierarchical amova showed that regions were important for explaining mtDNA variance in the Indian Ocean, but not Pacific, with social group sampling in the Atlantic too limited to include in analyses. Social groups were important in partitioning mtDNA and microsatellite variance within both oceans. Therefore, both geographic philopatry and social philopatry influence genetic structure in the sperm whale, but their relative importance differs by sex and ocean, reflecting breeding behaviour, geographic features and perhaps a more recent origin of sperm whales in the Pacific. By investigating the interplay of evolutionary forces operating at different temporal and geographic scales, we show that sperm whales are perhaps a unique example of a worldwide population expansion followed by rapid assortment due to female social organization.

Structures of K42N and K42Y sperm whale myoglobins point to an inhibitory role of distal water in peroxidase activity.

Sperm whale myoglobin (Mb) functions as an oxygen-storage protein, but in the ferric state it possesses a weak peroxidase activity which enables it to carry out H2O2-dependent dehalogenation reactions. Hemoglobin/dehaloperoxidase from Amphitrite ornata (DHP) is a dual-function protein represented by two isoproteins DHP A and DHP B; its peroxidase activity is at least ten times stronger than that of Mb and plays a physiological role. The `DHP A-like' K42Y Mb mutant (K42Y) and the `DHP B-like' K42N mutant (K42N) were engineered in sperm whale Mb to mimic the extended heme environments of DHP A and DHP B, respectively. The peroxidase reaction rates increased ∼3.5-fold and ∼5.5-fold in K42Y and K42N versus Mb, respectively. The crystal structures of the K42Y and K42N mutants revealed that the substitutions at position 42 slightly elongate not only the distances between the distal His55 and the heme iron but also the hydrogen-bonding distances between His55 and the Fe-coordinated water. The enhanced peroxidase activity of K42Y and K42N thus might be attributed in part to the weaker binding of the axial water molecule that competes with hydrogen peroxide for the binding site at the heme in the ferric state. This is likely to be the mechanism by which the relationship `longer distal histidine to Fe distance - better peroxidase activity', which was previously proposed for heme proteins by Matsui et al. (1999) (J. Biol. Chem. 274, 2838-2844), works. Furthermore, positive cooperativity in K42N was observed when its dehaloperoxidase activity was measured as a function of the concentration of the substrate trichlorophenol. This serendipitously engineered cooperativity was rationalized by K42N dimerization through the formation of a dityrosine bond induced by excess H2O2.

Low diversity in the mitogenome of sperm whales revealed by next-generation sequencing.

Large population sizes and global distributions generally associate with high mitochondrial DNA control region (CR) diversity. The sperm whale (Physeter macrocephalus) is an exception, showing low CR diversity relative to other cetaceans; however, diversity levels throughout the remainder of the sperm whale mitogenome are unknown. We sequenced 20 mitogenomes from 17 sperm whales representative of worldwide diversity using Next Generation Sequencing (NGS) technologies (Illumina GAIIx, Roche 454 GS Junior). Resequencing of three individuals with both NGS platforms and partial Sanger sequencing showed low discrepancy rates (454-Illumina: 0.0071%; Sanger-Illumina: 0.0034%; and Sanger-454: 0.0023%) confirming suitability of both NGS platforms for investigating low mitogenomic diversity. Using the 17 sperm whale mitogenomes in a phylogenetic reconstruction with 41 other species, including 11 new dolphin mitogenomes, we tested two hypotheses for the low CR diversity. First, the hypothesis that CR-specific constraints have reduced diversity solely in the CR was rejected as diversity was low throughout the mitogenome, not just in the CR (overall diversity π = 0.096%; protein-coding 3rd codon = 0.22%; CR = 0.35%), and CR phylogenetic signal was congruent with protein-coding regions. Second, the hypothesis that slow substitution rates reduced diversity throughout the sperm whale mitogenome was rejected as sperm whales had significantly higher rates of CR evolution and no evidence of slow coding region evolution relative to other cetaceans. The estimated time to most recent common ancestor for sperm whale mitogenomes was 72,800 to 137,400 years ago (95% highest probability density interval), consistent with previous hypotheses of a bottleneck or selective sweep as likely causes of low mitogenome diversity.

Can genetic differences explain vocal dialect variation in sperm whales, Physeter macrocephalus?

Sperm whale social groups can be assigned to vocal clans based on their production of codas, short stereotyped patterns of clicks. It is currently unclear whether genetic variation could account for these behavioural differences. We studied mitochondrial DNA (mtDNA) variation among sympatric vocal clans in the Pacific Ocean, using sequences extracted from sloughed skin samples. We sampled 194 individuals from 30 social groups belonging to one of three vocal clans. As in previous studies of sperm whales, mtDNA control region diversity was low (π = 0.003), with just 14 haplotypes present in our sample. Both hierarchical AMOVAs and partial Mantel tests showed that vocal clan was a more important factor in matrilineal population genetic structure than geography, even though our sampling spanned thousands of kilometres. The variance component attributed to vocal dialects (7.7%) was an order of magnitude higher than those previously reported in birds, while the variance component attributed to geographic area was negligible. Despite this, the two most common haplotypes were present in significant quantities in each clan, meaning that variation in the control region cannot account for behavioural variation between clans, and instead parallels the situation in humans where parent-offspring transmission of language variation has resulted in correlations with neutral genes. Our results also raise questions for the management of sperm whale populations, which has traditionally been based on dividing populations into geographic 'stocks', suggesting that culturally-defined vocal clans may be more appropriate management units.

Sperm whale population structure in the eastern and central North Pacific inferred by the use of single-nucleotide polymorphisms, microsatellites and mitochondrial DNA.

We use mitochondrial DNA (mtDNA) (400 bp), six microsatellites and 36 single-nucleotide polymorphisms (SNPs), 20 of which were linked, to investigate population structure of sperm whales (Physeter macrocephalus) in the eastern and central North Pacific. SNP markers, reproducible across technologies and laboratories, are ideal for long-term studies of globally distributed species such as sperm whales, a species of conservation concern because of both historical and contemporary impacts. We estimate genetic differentiation among three strata in the temperate to tropical waters where females are found: California Current, Hawai`i and the eastern tropical Pacific. We then consider how males on sub-Arctic foraging grounds assign to these strata. The California Current stratum was differentiated from both the other strata (P < 0.05) for mtDNA, microsatellites and SNPs, suggesting that the region supports a demographically independent population and providing the first indication that males may exhibit reproductive philopatry. Comparisons between the Hawai`i stratum and the eastern tropical Pacific stratum are not conclusive at this time. Comparisons with Alaska males were statistically significant, or nearly so, from all three strata and individuals showed mixed assignment to, and few exclusions from, the three potential source strata, suggesting widespread origin of males on sub-Arctic feeding grounds. We show that SNPs have sufficient power to detect population structure even when genetic differentiation is low. There is a need for better analytical methods for SNPs, especially when linked SNPs are used, but SNPs appear to be a valuable marker for long-term studies of globally dispersed and highly mobile species.

Alkylamine-ligated H93G myoglobin cavity mutant: a model system for endogenous lysine and terminal amine ligation in heme proteins such as nitrite reductase and cytochrome f.

His93Gly sperm whale myoglobin (H93G Mb) has the proximal histidine ligand removed to create a cavity for exogenous ligand binding, providing a remarkably versatile template for the preparation of model heme complexes. The investigation of model heme adducts is an important way to probe the relationship between coordination structure and catalytic function in heme enzymes. In this study, we have successfully generated and spectroscopically characterized the H93G Mb cavity mutant ligated with less common alkylamine ligands (models for Lys or the amine group of N-terminal amino acids) in numerous heme iron states. All complexes have been characterized by electronic absorption and magnetic circular dichroism spectroscopy in comparison with data for parallel imidazole-ligated H93G heme iron moieties. This is the first systematic spectral study of models for alkylamine- or terminal amine-ligated heme centers in proteins. High-spin mono- and low-spin bis-amine-ligated ferrous and ferric H93G Mb adducts have been prepared together with mixed-ligand ferric heme complexes with alkylamine trans to nitrite or imidazole as heme coordination models for cytochrome c nitrite reductase or cytochrome f, respectively. Six-coordinate ferrous H93G Mb derivatives with CO, NO, and O(2) trans to the alkylamine have also been successfully formed, the latter for the first time. Finally, a novel high-valent ferryl species has been generated. The data in this study represent the first thorough investigation of the spectroscopic properties of alkylamine-ligated heme iron systems as models for naturally occurring heme proteins ligated by Lys or terminal amines.

Quantifying the impact of dependent evolution among sites in phylogenetic inference.

Nearly all commonly used methods of phylogenetic inference assume that characters in an alignment evolve independently of one another. This assumption is attractive for simplicity and computational tractability but is not biologically reasonable for RNAs and proteins that have secondary and tertiary structures. Here, we simulate RNA and protein-coding DNA sequence data under a general model of dependence in order to assess the robustness of traditional methods of phylogenetic inference to violation of the assumption of independence among sites. We find that the accuracy of independence-assuming methods is reduced by the dependence among sites; for proteins this reduction is relatively mild, but for RNA this reduction may be substantial. We introduce the concept of effective sequence length and its utility for considering information content in phylogenetics.

Female philopatry in coastal basins and male dispersion across the North Atlantic in a highly mobile marine species, the sperm whale (Physeter macrocephalus).

The mechanisms that determine population structure in highly mobile marine species are poorly understood, but useful towards understanding the evolution of diversity, and essential for effective conservation and management. In this study, we compare putative sperm whale populations located in the Gulf of Mexico, western North Atlantic, Mediterranean Sea and North Sea using mtDNA control region sequence data and 16 polymorphic microsatellite loci. The Gulf of Mexico, western North Atlantic and North Sea populations each possessed similar low levels of haplotype and nucleotide diversity at the mtDNA locus, while the Mediterranean Sea population showed no detectable mtDNA diversity. Mitochondrial DNA results showed significant differentiation between all populations, while microsatellites showed significant differentiation only for comparisons with the Mediterranean Sea, and at a much lower level than seen for mtDNA. Samples from either side of the North Atlantic in coastal waters showed no differentiation for mtDNA, while North Atlantic samples from just outside the Gulf of Mexico (the western North Atlantic sample) were highly differentiated from samples within the Gulf at this locus. Our analyses indicate a previously unknown fidelity of females to coastal basins either side of the North Atlantic, and suggest the movement of males among these populations for breeding.

pH-induced equilibrium unfolding of apomyoglobin: substitutions at conserved Trp14 and Met131 and non-conserved Val17 positions.

A number of residues in globins family are well conserved but are not directly involved in the primary oxygen-carrying function of these proteins. A possible role for these conserved, non-functional residues has been suggested in promoting a rapid and correct folding process to the native tertiary structure. To test this hypothesis, we have studied pH-induced equilibrium unfolding of mutant apomyoglobins with substitutions of the conserved residues Trp14 and Met131, which are not involved in the function of myoglobin, by various amino acids. This allowed estimating their impact on the stability of various conformational states of the proteins and selecting conditions for a folding kinetics study. The results obtained from circular dichroism, tryptophan fluorescence, and differential scanning microcalorimetry for these mutant proteins were compared with those for the wild type protein and for a mutant with the non-conserved Val17 substituted by Ala. In the native folded state, all of the mutant apoproteins have a compact globular structure, but are destabilized in comparison to the wild type protein. The pH-induced denaturation of the mutant proteins occurs through the formation of a molten globule-like intermediate similar to that of the wild type protein. Thermodynamic parameters for all of the proteins were calculated using the three state model. Stability of equilibrium intermediates at pH ~4.0 was shown to be slightly affected by the mutations. Thus, all of the above substitutions influence the stability of the native state of these proteins. The cooperativity of conformational transitions and the exposed to solvent protein surface were also changed, but not for the substitution at Val17.

Time-dependent atomic coordinates for the dissociation of carbon monoxide from myoglobin.

Picosecond time-resolved crystallography was used to follow the dissociation of carbon monoxide from the heme pocket of a mutant sperm whale myoglobin and the resultant conformational changes. Electron-density maps have previously been created at various time points and used to describe amino-acid side-chain and carbon monoxide movements. In this work, difference refinement was employed to generate atomic coordinates at each time point in order to create a more explicit quantitative representation of the photo-dissociation process. After photolysis the carbon monoxide moves to a docking site, causing rearrangements in the heme-pocket residues, the coordinate changes of which can be plotted as a function of time. These include rotations of the heme-pocket phenylalanine concomitant with movement of the distal histidine toward the solvent, potentially allowing carbon monoxide movement in and out of the protein and proximal displacement of the heme iron. The degree of relaxation toward the intermediate and deoxy states was probed by analysis of the coordinate movements in the time-resolved models, revealing a non-linear progression toward the unbound state with coordinate movements that begin in the heme-pocket area and then propagate throughout the rest of the protein.